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We consider the dynamics of a movable mirror (cantilever) of a cavity coupled through radiation 
pressure to the light scattered from ultracold atoms in an optical lattice. Scattering from different 
atomic quantum states creates different quantum states of the scattered light, which can be distin- 
guished by measurements of the displacement spectrum of the cantilever. We show that for large 
pump intensities the steady state displacement of the cantilever shows bistable behaviour. Due to 
, atomic back-action, the displacement spectrum of the cantilever is modified and depends on the 

■ position of the condensate in the Brillouin zone. We further analyze the occurrence of splitting of 

' the normal mode into three modes due to mixing of the mechanical motion with the fluctuations 

^Nj , of the cavity field and the fluctuations of the condensate with finite atomic two-body interaction. 

The present system offers a novel scheme to coherently control ultracold atoms as well as cantilever 
dynamics. 
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INTRODUCTION 

<^ 

{^JQ' In recent years two distinct subjects, optical-micro cavities and nano-mechanical resonators have become entangled 
I ' . I experimentally by underlying mechanism of optical, radiation pressure forces. The coupling of mechanical and optical 
^ ' degrees of freedom via radiation pressure has been a subject of early research in the context of laser cooling [l|, 
Si S 13 and gravitational- wave detectors Q]- Recently there has been a great surge of interest in the application of 
radiation forces to manipulate the center-of-mass motion of mechanical oscillators covering a huge range of scales 
^ , from macroscopic mirrors in the Laser Interferometer Gravitational Wave Observatory (LIGO) project to nano- 
• mechanical cantilevers 0, SB [13, U, vibrating microtoroids 13, 14| membranes [l5j and Bose-Einstein condensates 
3 l3| ■ The quantum optical properties of a mirror coupled via radiation pressure to a cavity field show interesting 



I ■ similarities to an intracavity Kerr-like interaction [18j . Recently, in the context of classical investigations of nonlinear 



regimes, the dynamical instability of a driven cavity having a movable mirror has been investigated [19|. Theoretical 
work has proposed to use the radiation-pressure coupling for quantum non-demolition measurements of the light field 

In the field of quantum degenerate gases, standard methods to observe quantuni properties of ultracold atoms are 
based on destructive matter- wave interference between atoms released from traps [21[ . Recently, a new approach was 
^ • proposed which is based on all optical measurements that conserve the number of atoms. It was shown that atomic 
' quantum statistics can be mapped on transmission spectra of high-Q cavities, where atoms create a quantum refractive 
index. This was shown to be useful for studying phase transitions between Mott insulator and superfluid states since 
various phases show qualitatively distinct spectra New possibilities for cavity opto-mechanics by combining the 
tools of cavity electrodynamics with those of ultracold gases is the motivation of the present work. Here we show 
that different quantum states of ultracold gases in optical lattice, confined in a cavity can be distinguished by the 
■ steady state displacement spectrum of the movable mirror (cantilever). The atomic and cantilever back-action shifts 
I the cavity resonance. The laser-pump is shown to coherently control the dynamics of the mirror. Changing the pump 
J* • intensity, one can switch between stable and bistable regimes. Due to coupling between the condensate wavefunction 
. 5^ [ and the cantilever, mediated by the cavity photons, the cantilever displacement spectrum is continuously modified as 
I the condensate moves across the Brillouin zone. We also show that in the presence of atom-atom interactions, the 
5h ' coupling of the mechanical oscillator, the cavity field fluctuations and the condensate fluctuations (Bogoliubov mode) 
I leads to the splitting of the normal mode into three modes (Normal Mode Splitting) . 
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FIG. 1: Optomechanical realization of parametric coupling of a mechanical oscillator to a cavity optical mode and Bose- Einstein 
condensate trapped in an optical lattice inside the cavity. 



CANTILEVER DISPLACEMENT SPECTRA AS A PROBE OF QUANTUM PHASES OF ULTRACOLD 

ATOMS 

We consider an elongated cigar shaped Bose-Einstein condensate (BEC) of N two-level ^'^Rb atoms in the \F — 1 > 
state with mass m and frequency uja of the \F = 1 >— > \F' — 2 > transition of the D2 hne of ^''Rb, strongly 
interacting with a quantized single standing wave cavity mode of frequency tOc (Fig.l). The cavity field is also coupled 
to external fields incident from one side mirror. In order to create an elongated BEC, the frequency of the harmonic 
trap along the transverse direction should be much larger than one in the axial (along the direction of the optical 
lattice) direction. The system is also coherently driven by a laser field with frequency ujp through the cavity mirror 
with amplitude 77. The cavity mode is coupled to a mechanical oscillator (movable mirror) with frequency Qm via the 

dimensionless parameter e = {xo/^lm)-^\x=o- Here xq = y/h/2mef f^m is the zero point motion of the mechanical 
mode, nieff is its effective mass. It is well known that high-Q optical cavities can significantly isolate the system from 
its environment, thus strongly reducing decoherence and ensuring that the light field remains quantum- mechanical 
for the duration of the experiment. The harmonic confinement along the directions perpendicular to the optical 
lattice is taken to be large so that the system effectively reduces to one-dimension. This system is modeled by the 
opto-mechanical Hamiltonian (i?om) in a rotating wave and dipole approximation. 



where = cOp — LOa and Ac = ujp — ujc are the large atom-pump and cavity-pump detuning, respectively. Here 
cr+,(T~ are the Pauli matrices. The atom-field coupling is written as g{x) = gocos{kx). Here a and dm are the 
annihilation operators for a cavity photon and the mechanical mode respectively. The input laser field populates the 
intracavity mode which couples to the cantilever through the radiation pressure and the atoms through the dipole 
interaction. The field in turn is modified by the back-action of the atoms and cantilever. It is important to notice 
the nonlinearity in Eqn. ([T]) arising from the coupling between the intracavity intensity and the position quadrature 
of the cantilever. The system we are considering is intrinsically open as the cavity field is damped by the photon- 
leakage through the massive coupling mirror and the cantilever is connected to a bath at finite temperature. In 
the absence of the radiation-pressure coupling, the cantilever would undergo a pure Brownian motion driven by its 
contact with the thermal environment. Here we also assume that there is no direct coupling between the atoms 
and the cantilever, though this coupling could give rise to some interesting physics. Since the detuning Aa is large, 
spontaneous emission is negligible and we can adiabatically eliminate the excited state using the Heisenberg equation 

of motion = — [Hgrm o'^]- This yields the single particle Hamiltonian 




— ihg{x) [(T^a — a a 



HAcO,' a + hTlmOm drn 

] — iri{a — a)) + heflma'' a{am + aJ,J 



(1) 



Ho = -^ hAca^a + cos^{kx) [Vci{r) + hUoa''a] 

— iri{a — a^) + hflmam^ Om + heflma^ d{a„i + aj„). (2) 

g-2 

The parameter Uq — is the optical lattice barrier height per photon and represents the atomic backaction on 

the field Vci{r) is the external classical potential. Here we will always take Uq > 0. In this case the condensate 
is attracted to the nodes of the light field and hence the lowest bound state is localized at these positions which leads 
to a reduced coupling of the condensate to the cavity compared to that for Uq < 0. Along x, the cavity field forms 
an optical lattice potential of period A/2 and depth {hUo < a^a > +Vci). We now write the Hamiltonian in a second 
quantized form including the two body interaction term. 



H = J d^x^'<{r)Ho^{r) 



^147ra^ r ^3^^t(,^)^t(j;')^(^^(jr) (3) 
2 m J 

where \I'(r) is the field operator for the atoms. Here as is the two body s-wave scattering length. The corresponding 
opto-mechanical-Bose-Hubbard (OMBH) Hamiltonian can be derived by writing ^'(r) = bjw{r — rj), where w{f— 

rj) is the Wannier function and bj is the corresponding annihilation operator for the bosonic atom at the j*'* site. 
Retaining only the lowest band with nearest neighbor interaction, we have 

H^EoY, bib, +eY^ {b^.b, + 

j 3 

+ {hUoa^a + V,i) \ Jo % + jJ2 (^Vi^^ + [ + I ^ 

— hAcO^a — ihri{d — a^) + hflmOm^ Om + heVtmo) d{am + aj„) (4) 



where 



J.^[d^x.ir-r,).osHkZir-r,) 

The OMBH Hamiltonian derived above is valid only for weak atom- field nonlinearity ,24]. The nearest neighbor 
nonlinear interaction terms are usually very small compared to the onsite interaction and are neglected as usual. We 
now write down the Heisenberg-Langevin equation of motion for the bosonic field operator bj , the internal cavity 
mode d and the mechanical mode CI TIT QiS) 



bj = —i{Uad''d + 



h 
lEor 



) { JoS, + J(6j+i + - ^ + 6,_i} 



(6) 
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FIG. 2: Steady state displacement spectra of the cantilever. The single Lorentzian (thick line) for the MI reflects the non- 
fluctuating atom number. Many Lorentzians for SF (thin line) reflect atom number fluctuations. Here we have used gm/So = 0.2, 
^m/So = 30, Fm/So = O.OOf and k/Sq = 0.1 



(7) 



(8) 



Here n and Tm characterizes the dissipation of the optical and mechanical degree of freedom resectively. Here, 
we follow a semi-classical theory by considering noncommuting noise operators for the input field, i.e., (Cp(i)) = 0, 
{^li't')^p{t)) = "npSif - t), (fp(t')C|(0) = ("p + - t)i and a classical thermal noise input for the mechanical 

oscillator, i.e. {U{t)) = 0, {^Ut')U{t)) = {U{t')^L{t)) = n„,S{t' - t), in Eqs. 1^. The quantities and Up 
are the equilibrium occupation numbers for the mechanical and optical oscillators, respectively. We consider a deep 
lattice formed by a strong classical potential Vc\{v)^ so that the overlap between Wannier functions is small. Thus, we 
can neglect the contribution of tunneling by putting E = Q and J = . Under this approximation, the matter-wave 
dynamics is not essential for light scattering. In experiments, such a situation can be realized because the time scale 
of light measurements can be much faster than the time scale of atomic tunneling. One of the well-known advantages 
of the optical lattices is their extremely high tunability. Thus, tuning the lattice potential, tunneling can be made 
very slow [25| . 

The steady state value of the position quadrature Xm.s = o.m,s + a[,i s (the subscript s denotes the steady state 
value) is found as 



r2 



(9) 



where, 



(Ac - UqJqN - e^,nXrn,s? + ^2/4 

Here N = Vjhj. The equation for a\as is very important for the understanding the physics behind this problem. 
From it we clearly see how the cantilever and atom dynamics affects the steady state of the intracavity field. The 
coupling to the mirror and the atoms shifts the cavity resonance frequency and changes the field inside the cavity in a 
way to induce a new stationary intensity. The change occurs after a transient time depending on the response of the 



cavity and the strength of the couphng to the cantilever and the ultracold atoms. Following 23 Eqs. (|9ll0p allows 
to express Xm.s as a function /(ni, fiAf ) of atomic occupation number operators and calculate their expectation 
values for prescribed atomic states | ^t) . 

For the Mott state (fij) mi = Ij atoms are well localized at the jth site with no number fluctuations. It is represented 
by a product of Fock states, i.e. |5')mi = Yi^Li = kii Qm), with expectation values 



(/(ni, ...,nM))Mi = /(gi,...,(7M), (11) 
For simplicity we consider equal average densities {nj)Mi — N/M = no . 

In SF state, each atom is delocalized over all sites leading to local number fluctuations. It is represented 
by superposition of Fock states corresponding to all possible distributions of N atoms at M sites: |^)sf = 
J2qi qj^j /^yql\...qM^■\ql^ •■•7 1m)- Expectation values of light operators can be calculated from 



(/(ni,...,nAf))sF = 77]v Y] ' /(gi, gAf ), (12) 

^yiv ^ qi\...qM\ 
qi,---,qM 

representing a sum of all possible "classical" terms. Thus, all these distributions contribute to scattering from a 
SF, which is obviously different from (/(fii, ?t.a/))mi with only a single contributing term. 

In Fig. 2, we represent < \xm,s\ > as a function of Ac/(5o, where 5^ = UqJq. Clearly for the SF case, the displacement 
spectra is a sum of Lorentzians with different dispersion shifts. A comb like structure is seen if each Lorentzian is 
resolved. In the Mott state however, a single Lorentzian (thick line) is noticed. These were the typical structures 
predicted for the transmission spectra[23|. Since the coupling between the cantilever and the condensate wavefunction 
is mediated by the cavity field, we find that the displacement spectrum of the cantilever maps the distribution function 
of the ultracold atoms. An important condition to resolve the Lorentzians is k < Sq and this condition is easily met 
in present experiments [2^. Detecting the mirror's motion is straight-forward, since the optical phase shift is directly 
proportional to the mirror's displacement. Typically, the lorentizian frequency spectrum of the mirror's position is 
obtained in this way. The peak width yields the total damping rate, including the effective opto-mechanical damping. 
The area under the spectrum reveals the variance of the mirror displacement., which is a measure of the effective 
temperature. 



BISTABLE BEHAVIOUR 



We now consider the case of large number of atoms and hence treat the BEC within the mean field framework 
and assume the tight binding approximation where we replace bj by and look for solutions in the form of Bloch 
waves ipj = Ukexp{ikjd)exp{—ifit/h) [irt - Here /i is the chemical potential, d is the periodicity of the lattice and 

— ^ b\bj = no (atomic number density), M is the total number of lattice sites. From Eqns. ([9]) and pO|) . we obtain 
M f ■' 

a cubic equation in Xm,s- 

,3 _iA , (A^ + .V4) 8^^^ ^ 

Here, A = Ac — UoN{Jo + 2Jcos{kd)). Note that now wc do not ignore the tunneling term J. A plot of |a;m,s| 
versus A./60 for two different values of the pump parameter 77 is shown in Fig. 3. Clearly for higher pump intensity the 
system shows a bistable behaviour. For pump rates higher than a critical value, we find three steady state solutions for 
the mirror displacement, with two of them being stable. The system prepared below resonance will follow the steady 
state branch until reaching the lower turning point, where a non-steady-state dynamics is excited. This dynamics is 
governed by the time scale of the mechanical motion of the mirror because the cavity damping is almost two orders 
of magnitude faster (Fm << k). Fig. 3(b) shows the steady state displacement spectra |a;,„_s| as a function of the 
quasi-momentum for two different values of the cavity detuning Ac/Sq = 0.5 (thin line) and Ac/So — —0.5 (thick line). 
The cantilever phonons develop a quasi-momentum dependence due to strong coupling with the condensate mediated 
by the cavity photons. The atomic back-action modifies both the cavity field and the cantilever displacement. As 
the condensate moves across the Brillouin zone, the atom-field interaction changes and as a result the cantilever 
displacement spectra is continuously modified. For Ac/i5o = —0.5, a bistable behaviour is seen when the condensate 
is at the edge of the Brillouin zone (here A = 0). This bistable behaviour is absent when Ac/Sq — 0.5. The position 
of the condensate in the Brillouin zone is easily manipulated by accelerating the condensate. 



kd 



FIG. 3: (a): Steady state displacement spectrum as a function of A/5o for fc = 0. The thin curve corresponds to pump intensity 
77/(5o = 0.2 and the thick line corresponds to rj/5o = 1.0. For the larger pump intensity, we find three steady-state solutions, 
with two of them being stable. The system prepared below the resonance will follow the steady state branch until reaching 
the lower turning point, where a non-steady state dynamics is excited. The parameters taken are, eQ,m/5o ~ 0.2, flm/So — 30, 
Fm/^o ~ 0.001 and h/Sq = 0.08. (b): The displacement spectrum modified by the atomic backaction through the cavity 
photons as a function of the quasi-momentum (kd) for two different values of Ac/ So- For Ac/ So = —0.5, a bistable behaviour 
is seen (thick line) when the condensate is at the edge of the Brillouin zone (fed = ±7r and A = 0) . This bistable behaviour is 
absent when Ac/So — 0.5 (thin line). 

DYNAMICS OF SMALL FLUCTUATIONS: NORMAL MODE SPLITTING 

Here we show that the coupling of the mechanical oscillator, the cavity field fluctuations and the condensate fluctu- 
ations (Bogoliubov mode) leads to the splitting of the normal mode into three modes (Normal Mode Splitting(NMS)). 
The optomechanical NMS however involves driving three parametrically coupled nondegencrate modes out of equilib- 
rium. The NMS does not appear in the steady state spectra but rather manifests itself in the fluctuation spectra of 
the mirror displacement. To this end, wc shift the canonical variables to their steady-state values (i.e. d dg + d , 

dm dm,s + dm, bj — > (v^ + ^)) ^'^^ linearize to obtain the following Heisenberg-Langevin equations neglecting 

V M 

atomic losses due to heating: 

b= {ly + 2Uef f}b - iUef fb'' - igdd + d^) (14) 
d = (^iAd - ^) « - ("m + "L) - i9c{b + b^) + V^^pit), (15) 

dm = (^-i^rn - Sm ~ «^ (S -|- d''') + \/l\i^m{t). (16) 

Here, Ueff = 9c = UoJoVN\ds\ , v = C/q Jold^p + -^^^^ -f = - V^N - eQ^Xm^s is the detuning 

with respect to the renormalized resonance and A^ < leads to cooling. Note that Xm,s is negative and hence makes 
Ad positive. The atomic and cantilever back-action modifles the cavity detuning. The optomechanical coupling rate is 
given by gm = ^dge^m and |dsp gives the mean resonator occupation number. In deriving the above equation, we have 
assumed that ds to be real. As before, we assume negligible tunneling [J = E = 0) and hence we drop the site index 
j from the atomic operators.. We will always assume ^ k. Equations (|14ll5ll6p and their Hermitian conjugates 
constitute a system of six first-order coupled operator equations, for which the Routh-Hurwitz criterion implies that 



the system is stable only for g„i < y i^dA'^ + K^/4)f7„i/| A^l and AdA^ > 0, where A^ = Ad — 2gc/{v + SUeff). 

We transform to the quadratures: Xm = dm + dm, Pm = iidm ~ Om), Xp = d + d^, Pp = i(d^ — d), Xh = b + b\ 
Pb = i{w —b). The displacement spectrum in Fourier space for rij, = is found as 
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FIG. 4: (a) : Normalized plot of the displacement spectrum Sx (io) for two values of the atomic two-body interaction, Uef / /f^m = 
(thin line) and Ueff/^m ~ 0.3 (thick line). Parameters chosen are: Qm/^m ~ 0.4, n/Qm = 0.1, rm/^m = 0.01, Qc/^m = 0.1, 
Ad/r^m = —1 and v/Q.m ~ 0.01. A clear NMS is observed in the presence of atom-atom interaction. , (b): Normalized plot of 
the displacement spectrum Sx{i^) for the same parameters as (a) but with a stronger atom-photon coupling gc/0,m = 0.2. The 
NMS is now much more prominent compared to that in (a). 
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where, 



^2/4 



r.(^) 



(17) 



(18) 



fis(t^) = 



and 



(19) 



r.(c^) = - 



{T8 + K^uj^nt} 



(20) 



(- + A2 - ^2)(^2 „ ^ _^ 3f^^^^.)) _ 4^^^2(^ ^ 



(21) 



(22) 



This spectrum is characterized by a mechanical susceptibility xi'-^) that is driven by thermal noise (cx rim) and by the 
quantum fluctuations of the radiation pressure (quantum back-action) and quantum fluctuations of the condensate. 
Note that in the absence of the atom-photon couphng {gc = 0), the displacement spectrum EqnfTTl reduces to that 
found in 

In Fig. 4a, we show a normalized plot of the displacement spectrum Sx{io) for two values of the atomic two-body 
interaction, Ueff/^m — (thin line) and Ueff/^m — 0.3 (thick line). Parameters chosen are: gm/^m = 0.4,K/r2„i = 
0.1,r„j/i^m = O.OI, gd^m ~ 0.1, Ad/fim = ^ 1 and v/Vlm = 0.01. In the absence of the interactions, we observe the 
usual normal mode spliiting into two modes and we flnd that due to a flnite atom-atom interaction the normal mode 
splits up into three modes. The NMS is associated with the mixing between the mechanical mode and the fluctuation 
of the cavity field around the steady state and the fluctuations of the condensate (Bogoliubov mode) around the mean 
field. The origin of the fluctuations of the cavity field is the beat of the pump photons with the photons scattered 



from the condensate atoms. The frequency of the Bogohubov mode in the low momentum hmit is w ^Ueff- Hence 
in the absence of interactions, the Bogohubov mode is absent and the system simply reduces to the case of two mode 
coupling (i.e coupling between the mechanical mode and the photon fluctuations). In the presence of finite atom- 
atom interaction, the mechanical mode, the photon mode and the Bogohubov mode forms a system of three coupled 
oscillators. It is important to note that splitting in the displacement spectrum occurs only for g„i > due to 

finite width of the peaks. Fig. 4(b) illustrates the influence of increasing the atom-photon coupling gc — 0.2. The 
NMS is now much more prominent compared to the case with gc ~ 0.1. Similar three coupled oscillator experimental 
results where two coupled cavities, each containing three identical quantum wells [29] and one micro-cavity containing 
two quantum wells ^3Q] have been reported. One experimental limitation could be spontaneous emission which leads 
to momentum diffusion and hence heating of the atomic sample 17| . The parameters chosen here are within current 
technological reach 1^, 17, 31| and hence the observation of the paramertic NMS is within experimental reach. 



CONCLUSIONS 



In summary we have analyzed a novel scheme of cavity-opto-mechanics with ultracold atoms. We showed that the 
steady state displacement spectra of a cantilever coupled indirectly to a gas of ultracold atoms in an optical lattice 
through the cavity field are distinct for different quantum phases of equal densities. Further, we showed that in the 
mean-field of the atoms, the cantilever displacement shows a bistable behviour for high pump intensities and exhibits 
a dependence on the position of the condensate in the Brillouin zone due to the back-action of the condensate on the 
cantilever. The cantilever displacement spectrum shows a bistable behaviour when the condensate is at the edge of 
the Brillouin zone. In the presence of atom-atom interactions, the coupling of the mechanical oscillator, the cavity 
field fluctuations and the condensate fluctuations (Bogohubov mode) leads to the splitting of the normal mode into 
three modes (Normal Mode Splitting). The system described here shows a complex interplay between distinctly three 
systems namely, the nano-mechanical cantilever, optical micro-cavity and the gas of ultracold atoms. The quantum 
state of the degenerate gas influences the cantilever displacement spectra via the cavity photons, while on the other 
hand the cantilever displacement modifies the cavity field which in turn modifies the properties of the condensate. 
Manipulating the position of the condensate in the Brillouin zone and the atomic tunneling, we now have a new 
possibility to efficiently control the cantilever dynamics and at the same time the cantilever dynamics can be used to 
control the atomic dynamics. 
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